A = Point(yAxis)
B = Point(xAxis)
b = Line(A, B)
C = Point(zAxis)
D = Point(xAxis)
a = Line(C, D)
h = Line(B, a)
a_1 = Plane(b, h)
i = OrthogonalLine(C, a_1)
j = OrthogonalLine(D, a_1)
b_1 = Plane(i, j)
k = IntersectPath(b_1, a_1)
P = Intersect(k, b)
l = OrthogonalLine(P, a, space)
Q = Intersect(l, a)
α = Angle(l, a)
β = Angle(l, b)